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ABSTRACT 

Context. Supernova 1987 A (SN 1987 A) in the neighbouring Large Magellanic Cloud offers a superb opportunity to follow the evolu- 
tion of a supernova and its remnant in unprecedented detail. Recently, far-infrared (far-lR) and sub-mm emission was detected from 
the direction of SN 1987 A, which was interpreted as due to the emission from dust, possibly freshly synthesized in the SN ejecta. 
Aims. To better constrain the location and hence origin of the far-lR and sub-mm emission in SN 1987A, we have attempted to resolve 
the object in that part of the electro-magnetic spectrum. 

Methods. We observed SN 1987 A during July-September 2011 with the Atacama Pathfinder Experiment (APEX), at a wavelength 
of 350 yum with the Submillimetre APEX Bolometer CAmera (SABOCA) and at 870 yum with the Large APEX BOlometer CAmera 
(LABOCA). The 350-yum image has superior angular resolution (8") over that of the Herschel Space Observatory 350-/jm image 
(25"). The 870-yum observation (at 20" resolution) is a repetition of a similar observation made in 2007. 

Results. In both images, at 350 and 870 /im, emission is detected from SN 1987A, and the source is unresolved. The flux densities in 
the new (2011) measurements are consistent with those measured before with Herschel at 350 //m (in 2010) and with APEX at 870 
//m (in 2007). A higher dust temperature (a 33 K) and lower dust mass might be possible than what was previously thought. 
Conclusions. The new measurements, at the highest angular resolution achieved so far at far-lR and sub-mm wavelengths, strengthen 
the constraints on the location of the emission, which is thought to be close to the site of SN 1987A and its circumstellar ring structures. 
These measurements set the stage for upcoming observations at even higher angular resolution with the Atacama Large Millimeter 
Array (ALMA). 
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1. Introduction 

' SN 1987A offers us the unique opportunity to follow the early 
, evolution of a supernova (SN) and its transition into a super- 
nova remnant (SNR) in exquisite detail. Located in the Large 
Magellanic Cloud (LMC) at a relatively accurately known dis- 
. tance of 50 kpc, it is seen in a direction that suffers from little 
' interstellar extinction. While SN 1987 A was classified as type II 
" (Arnett et al. 1989), it showed distinct propeities with respect 
to typical plateau events, most probably linked to its unexpected 
■ progenitor star: the blue supergiant Sk -69° 202. 

SN 1987A's emerging remnant is being shaped by the im- 
pact of the initial blastwave and the expanding hot plasma upon 
an equatorial ring of matter, 1 "6 in diameter, deposited probably 
by the supernova progenitor star when it was still a red super- 
giant (Bouchet al. 2006; Dwek et al. 2010; Larsson et al. 201 1). 
The ejecta themselves are due to reach the ring in the near fu- 
ture. The increasing strength of the shock is also reflected in an 
exponential increase in radio brightness (Zanardo et al. 2010). 

Small amounts (Mj < 10""' Mq) of warm dust formed in 
about two years since the explosion (Suntzeff & Bouchet 1990; 
Bouchet, Danziger & Lucy 1991). This offers a great opportunity 
to study fresh supernova dust, because in older remnants it be- 
comes increasingly difficult to distinguish between dust formed 
in the ejecta or dust swept-up from the suiToundings (e.g., Rho 
et al. 2008, 2009; Sandstrom et al. 2009; van Loon et al. 2010; 
Otsuka et al. 2010; Barlow et al. 2010). Early detection also al- 
lows one to monitor the evolution of the dust, which might be al- 



tered or even destroyed before enteiing the interstellar medium. 
Claims of observaions of dust formation have been made in more 
distant SNe (e.g., Elmhamdi et al. 2003; Fox et al. 2009; Sakon 
et al. 2009; Inserra et al. 2011), though often SNe are seen to 
destroy pre-existing dust (e.g., Botticella et al. 2009; Wesson et 
al. 2010; Andrews et al. 201 1). SN 1987A offers prospects of lo- 
cating the cold dust directly, by spatially resolving its emission. 

Far-IR emission was detected from SN 1987 A with the 
Herschel Space Observatory by Matsuura et al. (201 1), who sug- 
gested a colder dust component (17-23 K) than that in the equa- 
torial ring (a: 170 K). They placed the cold dust within the ejecta, 
to reconcile the large mass of 0.5 Mq they deiive. This would 
imply enormous growth since the earlier detections. However, 
the Herschel images do not resolve the equatorial ring from the 
ejecta, and still allow for the dust to reside outside the ring - 
be it of progenitor origin or interstellar dust echoing the super- 
nova (Meikle et al. 201 1). We also remind the reader that chance 
superpositions do occur (e.g., van Loon & Oliveira 2003). 

SN 1987 A was observed at sub-mm wavelengths (870 fim) 
with the Large APEX BOlometer CAmera (LABOCA; Siringo 
et al. 2009) on the Atacama Pathfinder Experiment (APEX) 
for the first time in 2007 (Lakicevic et al. 2011). We observed 
this object again in 201 1 with LABOCA at 870 yum, and at 350 
fim with SABOCA (Submillimetre APEX Bolometer CAmera; 
Siringo et al. 2010), at three times higher angular resolution than 
the Herschel observations at that wavelength. We here present 
the results. These new observations place important constraints 
on the maximum spatial extent of the dust emission, and they 
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Fig. 1. APEX/LABOCA image of SN 1987A at 870 iim, beam 
FWHM 20"4. Contours are used where the brightness exceeds 
30 mJy beam-'. SN1987Ais located at 5^35™28?0,-69°16'll". 



provide an important zero-spacing value for planned observa- 
tions with the Atacama Large Millimeter Array (ALMA). 



2. Observations 

Observations of SN 1987A (programme ESO088.D-0252, P.I. 
Masa Lakicevic) were made with APEX, a 12-m single-dish sub- 
mm telescope at the ALMA site at an altitude of 5100 m. We 
used the bolometer array cameras LABOCA and SABOCA, in 
July and September 201 1 under excellent conditions - the water 
vapour column was < 0.2 mm for most of the SABOCA obser- 
vations and < 0.5 mm for the LABOCA observations. 



2.1. LABOCA observations at 870 ^im 

SN 1987 A was observed with LABOCA at a wavelength of 870 
urn (345 GHz) within a 60-GHz wide band (Av/v ^ 1/6). The 
on-source integration time was 5 h (6.1 h including overheads). 
The field was imaged in spiral raster map mode. The beam has a 
full-width at half-maximum (FWHM) of ~ 19.2". 

The data were reduced with the bolometer array analysis 
(BoA) software. The data were first corrected for atmospheric 
attenuation determined through "skydips", applying a correction 
derived from observations of planets and secondary flux calibra- 
tors taken at similar elevation as the target. Then several rounds 
of correlated flux variation ("skynoise") removal, despiking, re- 
moval of noisy bolometers, and baseline subtractions were ap- 
plied, before maps of the individual scans were created. These 
were then combined. The resulting map was then smoothed, and 
values below zero set to zero, and used as a source model for the 
next iteration. The model was subtracted from the data before 
skynoise, spike and baseline removal and added back before cre- 
ating the next map. In total 35 iterations were performed, with in- 
creasing skynoise, spike and baseline removal settings while re- 
ducing the smoothing. This procedure allows one to recover ex- 
tended structures, reduces the negative "bowls" seen frequently 
in bolometer on-the-fly maps, and reduces the noise, while yield- 
ing a final map at the full spatial resolution of the instrument. 

The image was smoothed with half a beam; the resulting 
beam size was 20!'4. It is presented in full in Fig. 1, and is 
zoomed-in on the position of SN 1987A in Fig. 2 (top right). 



The source appears to sit on filamentary emission extending to- 
wards the south-west, similar to that reported by Lakicevic et al. 
(2011). 

The flux density of the central point source is F^qo - 19.6 + 
5.6 mJy, where the uncertainty combines the rm.s. noise level of 
4 mJy (8 mJy in the unsmoothed map) and LABOCA's absolute 
calibration uncertainty of ~ 20%. This brightness is consistent 
within the errors with the result from Lakicevic et al. (20 11), viz. 
21+4 mJy, which was obtained from measurements in 2007. 



2.2. SABOCA observations at 350 /jm 

SABOCA is a continuum receiver that consists of 39 bolometers 
operating at a wavelength of 350 fim (857 GHz) within a 120- 
GHz wide band (Av/v ^ 1/7). The on-source integration time 
was 15.6 h (22.7 h including overheads). The field was imaged 
in spiral raster map mode. The beam has a FWHM of 7"8 + 0''5. 

The data reduction was performed with two different pack- 
ages - BoA and crush Comprehensive Reduction Utility for 
SHARC-2; Kovacs 2008 - which use different algorithms. 

In CRUSH, we used the faint source reduction option, with a 
source size filtering of 10", and ran 20 iterations. The image, 
smoothed with half a beam, is presented in Fig. 2 (top left). 
Since the pointing in September was slightly off'set from that in 
July, the data from July and September were reduced separately 
and then the September data were shifted to match the July data 
and averaged. The resulting beam has an eff'ective FWHM of 
8!'4. The flux density of the unresolved source at the position of 
SN 1987A on the combined image is F350, crush = 58 + 13 mJy, 
where the uncertainty combines the r.m.s. noise level of 6 mJy 
and SABOCA's absolute calibration uncertainty of ~ 20%. 

The data reduction strategy within BoA corresponds to the 
one used for LABOCA, except that the atmospheric opacities 
were determined from measurements of the amount of precip- 
itable water vapour (PWV) done in parallel to the observa- 
tions, using an atmospheric model (ATM). Because the cali- 
bration measurements were made at the same time as the ob- 
servations and at the same azimuth and elevation, they yield 
a more accurate opacity correction. The image obtained using 
BoA, smoothed with half a beam, is presented in Fig. 2 (top, cen- 
tre). Here, the July and September data were aligned and placed 
at the centre of the image before the reduction. The measured 
flux density of SN 1987 A is now F350, boA -44 + 1 mJy. 

The maps reduced with BoA and crush diff'er slightly in the 
flux level, probably due to the different calibration strategies 
(crush uses opacities derived from skydips). The flux densities 
are consistent within the calibration uncertainties of ~ 20-30%. 
The extended emission in the BoA reduced image is not evident 
in the crush reduced map, because here the algorithm is opti- 
mised to remove correlated signal variations from the data time 
streams, which includes extended emission structures. 

The images from the two reduction techniques both confirm 
that the object at the position of SN 1987 A is a point source. The 
detection is of much higher significance than the uncertainty in 
the flux density suggests, because the latter includes systematic 
uncertainties mainly in the flux calibration. For the flux density 
we adopt the average from the two reductions, i.e. ^350 = 51 + 10 
mJy, where the uncertainty is the average from the individual 
uncertainties as they are strongly correlated. This is remarkably 
close to the Herschel measurement, viz. 49.3 +6.5 mJy, in 2010. 
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Fig. 2. Top: APEX images of SN 1987A - left: at 350 /im reduced with crush, FWHM 8"4; centre: at 350 //m reduced with BoA, 
FWHM 8"4; right: at 870 /im, FWHM 20;'4. Bottom: previous images from Herschel (Matsuura et al. 201 1) and APEX (Lakicevic 
et al. 201 1) - left: at 100 yum, FWHM 6'.'%; centre: at 350 /^m, FWHM 24:'9; right: at 870 //m with APEX, FWHM 23". SN 1987A 
is located at 5'>35'"28?0, -69°16'11". 



3. Discussion 

Models predict ~ 0.1-1 M© of dust to form in the ejecta of core- 
collapse SNe (Nozawa et al. 2003). Warm dust (» 170 K) was 
identified with the equatorial ring around SN 1987 A by direct 
imaging (Bouchet et al. 2006). Dust formed within the ejecta 
might be colder The cooling rate of a spherical grain of diam- 
eter a and temperature Tj is its luminosity, Lj = An{al2)-crT'^, 
resulting in a heat loss -mCdTd/dt where m is the mass of the 
grain and C its specific heat. Integrating, noting that the initial 
temperature rapidly vanishes, leads to 



Td 



C\'" 



la pL\ 
ll8crj 



:-l/3 



(1) 



where p is the grain's density. For p ~ 3 g cm"^, C ~ 1 J K ' 
g ' and a = 0.1 /im one finds that the grain cools to the tem- 
perature of the cosmic microwave background (3 K) within a 
matter of hours. First, collisions with the gas prevent this from 
happening, but as the ejecta expand the gas cools and the gas and 
grains detach (when the density falls below ~ 10^ cm"^). The in- 
terstellar radiation field (ISRF) keeps the dust at a temperature 
10 < Td < 30 K. This led Matsuura et al. (201 1) to suggest that 
the ~ 20-K dust attributed to the fai--IR emission from SN 1987 A 
formed within the SN ejecta. 



3. 1 . Spectral energy distribution (SED) 

In Fig. 3 we present the SED of SN 1987A. Lakicevic et al. 
(2011) described the far-IR/sub-mm emission using a modified 
black body with a dust emission coefficient k/kq - (A/Ao)"; they 
found a = -1.5 and T^ x 18 K, similar to the values reported 
in Matsuura et al. (2011) and hence a similar dust mass (a; 0.5 
Mq) depending on kq. Synchrotron emission contributes ~ 5 and 
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Fig. 3. SED of SN 1987A. Data: Biermannet al. (1992);Zanardo 
et al. (2010); Matsuura et al. (2011); Lakicevic et al. (2011, 
2012); this work. 



~ 9 mJy at 0.35 and 0.87 mm, respectively, as estimated from 
extrapolation of a powerlaw Fy oc A" with slope a - -(-0.7 (Fig. 
3). Limits on the contribution from free-free emission at those 
wavelengths are rather more stringent (Lakicevic et al. 201 1). 
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The new measurements at 0.35 and 0.87 mm show no change 
from the previous measurements at those wavelengths with 
Herschel (Matsuura et al. 2011; data from August 2010) and 
APEX (Lakicevic et al. 201 1 ; data from July 2007), respectively. 
Hence there is no evidence for an evolution of the dust on those 
timescales, be it the amount, temperature, or grain sizes. 

However, values of a closer to zero may not be unexpected 
at sub-mm wavelengths (Paradis et al. 2011). Adopting a - Q, 
and raising the dust temperature to 33 K results in a better fit 
to the 0.1, 0.35- and 0.87-mm measurements (Fig. 3), leaving 
only the Herschel 160- and 250-/im measurements discrepant. A 
compromise is obtained for a - -0.5, at 27 K, which is at the 
extreme end of the models in Fig. 7 of Paradis et al. (201 1); this 
reconciles the 160-/im datum (Fig. 3). The Herschel data could 
be affected by line emission close to SN 1987 A: the 160-;um 
band includes the [Cii] 158-/im and [Nii] 205-;um lines, while 
the 250-/im band includes three highly excited '^CO lines with 
another '^CO line and the [N ii] line at the edges of the band. In 
warm gas, these '^CO transitions can be brighter than those at 
longer wavelengths, and the 350- and 500-/im bands each only 
include two '^CO lines, which additionally reduces the contri- 
bution from line emission at A> 300 /im. The fit for a = (and, 
though to a lesser degree, the fit for a = -0.5) agrees with the 
measurement at 1.3 mm made in 1990 (Biermann et al. 1992), 
hinting that the cold dust has been unaltered since day 1290 and 
that it might well have already been in place prior to the explo- 
sion. Following Evans et al. (2003) we obtain 



Md/Mo=1.5xlOVmr 



(2) 



where k^ cannot exceed the ratio of grain surface area and grain 
mass. Hence, a spherical black body with a = 0.1 /im and p = 3 
g cm"^ would yield /c = 2 x 10^ cm^ g"' . This would reduce the 
mass estimate to M^ ~ 0.008 M© and within the realms of the 
formation through mass loss from the SN progenitor (assuming 
several M© were lost and a gas-to-dust ratio of a few hundred). 
For a - -0.5 or -1.5 at /I > 100 /im we would have found 
Md ~ 0.024 or 0.24 Mq, respectively. Mennella et al. (1998) find 
K ~ 10-10^ cm^ g"' at far-IR-mm wavelengths, depending on 
grain composition, but smaller grains, needles, or fluflFy grains all 
increase the surface-to-mass ratio and hence k. Indeed, for SNe 
of type lib - probably a good match to SN 1987 A - Nozawa 
et al. (2010) predict grains of only 0.001-0.003 //m size, and 
Wickramasinghe & Wickramasinghe (1993) explained the mm 
radiation from SN 1987 A in 1990 by iron whiskers with k ~ \(f 
cm^ g"' at far-IR-mm wavelengths. 

3.2. Where is the dust? 

Our APEX 350-/im image of SN 1987 A has a beam size (area) 
a tenth of that of the Herschel 350-/zm image (Matsuura et al. 
201 1). The emission from SN 1987A remains unresolved, i.e. it 
is located well within 8" (corresponding to 2 pc at the distance of 
the LMC of 50 kpc) from the explosion site. Even if it were due 
to a light echo (Meikle et al. 201 1) this suggests the echoing dust 
is probably situated < 2 pc from the explosion site and was most 
likely produced by the progenitor of SN 1987 A. Little room now 
remains for cold dust to surround the equatorial ring, i.e. in the 
dense wind material from the progenitor; the dust is more likely 
to reside within the ring or inside of it (e.g., in the ejecta). 

The grains might be heated to 33 K by the decay of ^''Ti 
within the ejecta (Jerkstrand, Fransson & Kozma 2011). Small 
grains would naturally be hotter within the ISRF, and also more 
prone to stochastic heating, which could imply that not all dust 
mass is accounted for (the addition of a much colder component 



could mimic a flatter long-wavelength tail). Small grains might 
result from sputtering by the hot gas and X-rays impinging upon 
the ring or even the ejecta (Larsson et al. 201 1). While Dwek et 
al. (2010) did find an additional population of small grains in the 
equatorial ring, these were much hotter (~ 350 K). Hence the 
cold dust might be most easily explained as due to condensation 
within the ejecta, but there are scenarios possible that reduce the 
inferred mass to «c 0.5 Mq, leaving other options open. 

4. Conclusions 

We presented new images of SN 1987 A at wavelengths of 350 
and 870 fan, obtained with the APEX telescope, showing very 
clear detections of largely unresolved emission. The image at 
350 fim improves the constraints on the location of the emis- 
sion by an order of magnitude in area over that obtained from 
Herschel observations, to well within 8" from the explosion 
site. The far-IR and sub-mm SED is consistent with emission 
from cold dust, but the temperature may be as high as ~ 33 K 
if the emission coefficient is wavelength-independent - though 
this would require a contribution from line emission to the 
flux within the 160- and 250-/im Herschel bands. Smaller, non- 
spherical, or fluffy grains would reduce the inferred dust mass, 
possibly enough to become consistent with pre-SN production. 
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